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ABSTRACT 

A growing body of evidence suggests that multiple dynamo mechanisms can drive magnetic vari¬ 
ability on different timescales, not only in the Sun but also in other stars. Many solar activity proxies 
exhibit a quasi-biennial (~2 year) variation, which is superimposed upon the dominant 11 year cy¬ 
cle. A well-characterized stellar sample suggests at least two different relationships between rotation 
period and cycle period, with some stars exhibiting long and short cycles simultaneously. Within 
this sample, the solar cycle periods are typical of a more rapidly rotating star, implying that the 
Sun might be in a transitional state or that it has an unusual evolutionary history. In this work, we 
present new and archival observations of dual magnetic cycles in the young solar analog HD 30495, 
an ~1 Gyr-old G1.5V star with a rotation period near 11 days. This star falls squarely on the rela¬ 
tionships established by the broader stellar sample, with short-period variations at ~1.7 years and a 
long cycle of ~12 years. We measure three individual long-period cycles and find durations ranging 
from 9.6-15.5 years. We find the short-term variability to be intermittent, but present throughout 
the majority of the time series, though its occurrence and amplitude are uncorrelated with the longer 
cycle. These essentially solar-like variations occur in a Sun-like star with more rapid rotation, though 
surface differential rotation measurements leave open the possibility of a solar equivalence. 

Subject headings: stars: individual HD 30495 — stars: activity — stars:solar-type — stars: rotation 
— stars: dynamo 


1. BACKGROUND AND CHARACTERIZATION 

Stellar magnetic activity cycles have been known in the 
Sun since Schwabe and were shown to exist in other stars 


by Wilson (1978) using an activity index derived from 


flux measurements in the Ca II HK line cores performed 
at the Mount Wilson Observatory. The Mount Wilson 
survey from 1966-2003 remains the largest and longest 
campai gn inves tigating stellar activity, culminating with 


Baliunas et al. (19951 reporting on the variability of 111 
stars, 52 of whicK were found to demonstrate periodic 
behavior in a manner similar to our Sun. This periodic 
behavior is believed to be caused by a magnetic dynamo 
driven by rotational motions of plasma in the stellar inte¬ 
rior. Using the Mount Wil son data along with mea sure¬ 


ments of surface rotation, Saar & Baliunas (19921) and 


Soon et al. (1993) first reported two distinct brandies of 
cycling stars, the “active” and “inactive” branches dis¬ 
tinguished by their mean activity level and number of ro¬ 
tations per cycle, and furthermore found some stars with 
multiple prominent periodicities fit on both branches. 


the stellar sample of 

Saar & Brandenburg ( 

1999 

), 

Bohm- 
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Vitense (2007) hypothesized that these classes are due to 
separate dynamo mechanisms identifiable by their inte¬ 
rior shear layer, with the slow-rotating “inactive” branch 
operating near the base of the convective envelope, and 
the fast-rotating “active” branch working in rotational 
shear layers closer to the surface. Curiously, the Sun and 
its 11 yr cycle appear to be a unique outlier in this small 
stell ar sam ple, f alling between the two activity branches. 


Bohm-Vitense (2007) suggested that the Sun could be 
in transition from one dominant dynamo mechanism to 
another. 

In addition to the 11 yr solar cycle, short-term quasi- 
periodic variability has been observed in a number of so¬ 
lar phenomena. Various manifestations of the so-called 
q uasi-biennial oscillat i on (Q BO) of 0.6-4 yrs are reviewe d 


Bazilevskaya et al. (2014) and McIntosh et al. (2015). 


QBOs are found in records of sunspot number and area 
magnetic field measurements, solar irradiance, and in 
magnetically-sensitive phenomena such as filaments in 
H-alpha and variations of field-sensitive lines such as Ca 
II and Mn I. The QBO also appears in eruptive phe¬ 
nomena - flares, coronal mass ejections, and solar en¬ 
ergetic particle events - which arise from magnetically 
active regions. Shor t period variations hav e also been ob¬ 
served in the stars: Baliunas et al. (1995) reported nine 
stars with significant “secondary cycles” . Six of these 
(and one new addition) were part of the high-q uality 
activity cycles of the Saar & Brandenburg (1999) sam¬ 
ple, wit h the secondary cy cle falling on the “inactive” 


branch. Olah et al. (2009) performed a time-frequency 
analysis oiYnuiti-decadal photometry and Ca II emission 
for 20 stars and found 15 of them to exhibit multiple 
cycles. High-cadence SMARTS HK observations have 
found short-perio d variations (1.6 yr) on l Horologi (Met- 
calfe et aTpOlO) and e Eridani (2.95 yr), another dual- 
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cycle star w i th a long-term cyc le of 12.7 yr (Metcalfe 


et al. 2013). Fares et al. (20091 used Zeeman Doppler 


Imaging to observe a polarity-flipping cycle of ~2 years 
in the fast-rotating F6 star r Bootis, which had weak 
indications of a long -period activity cycle of 11.6 yr in 
Baliunas et al. (1995). This short-period cycle is distinct 
from the solar QBO phenomena, which does not reverse 
magnetic polarity. 

The origin of the solar QBO and its relation to the solar 
cycle is not understood. The discovery of periodic vari¬ 
ations of 1.3 yrs in the differential r otation of th e dee 
interior revealed by helioseismology (Howe et al. 2000 
suggests that the QBO is sub-surface m origin and may 
be an additional feature of a deep-interior dyn amo pro¬ 


cess respons ible for the 11-year cycle (see |Bazilevskaya 


et al. 


2014 and references therein). McIntosh et al. 


f2(J15) also points to a deep interior process, interring 
that this short-period variability is driven by the inter¬ 
action of two oppositely-signed magnetic activity bands 
deep in the interior of each hemisphere. Another pos¬ 
sibility is the distinct-dynamo scenario described above 
as an explanation for the two activity branches i n the 
Saar fe Bra ndenburgl fll999MBdh m-Vitense| ((2007| stel- 
iar sample. FlefcEeref^rll 2010j) find a ~2 yr variation 
in the frequency shift of solar p-niode oscillations, and 
locate the origin of the variations to be below the source 
of the 11-yr signal in the data. They hypothesize that 
spatially distinct dyn amo processes may be r esponsible 
for this phenomenon, 
a non-dynamo origin 


Chowdhury et al. (2009) suggested 
or QBOs: an instability caused by 
Rossby waves interacting with the tachocline. 

While stellar observations cannot match the level of de¬ 
tail in which the solar QBO is observed, the varied physi¬ 
cal conditions present in other stars may have an impact 
on the manifestations of these short-period oscillations 
that sheds additional light on their origins. In this work, 
we present new observations of variability in HD 30495 
(58 Eri), a nearby young solar analog that demonstrates 
both a long-term activity cycle and a short-period os¬ 
cillation, which may be analogous to the solar cycle and 
QBO. The upper section of Table [l] summarizes the mea¬ 
sured global properties of HD 30495. Photometry and 
spectroscopy show the star to be essentially solar-like, 
leading previou s authors to study this star as a poten- 


tial solar twin jCayrel de Strobel 1996 Porto de Mello 


et al.|[2014 ). Gaidos et al. 1 2000) found rotational mod¬ 


ulations in high-cadence Stromgren b and y photometry 
to determine a rotation period of 11.3 days, roughly 2.3 
times faster than the Sun. Due to the process of magnetic 
braking, older stars have slower rotations, giving t he well- 
known P rot ~ t ~ 0 5 age-rotation relationship (Skumanich 


1972). The faster rotation of HD 30495 thereby implies 


it is younger tha n the S un , and by the age-rotation rela¬ 
tionship given in Barnes|(2007) we obtain an age of about 
~1 Gyr. Observations of excess infrared flux attributed 
to a diffuse and distant debris disk of ~73 Earth-masses 


leftover from formation (Habing et al. 2001) give fur¬ 
ther evidence for a young age. Spectroscopic searches 
of similar “Vega-like” main-sequence objects with excess 
infrared emission have ruled out the possib ility of dense 


concentrat ions of gas close to the star (Habing et al. 2001 
Liseau||1999|). Based on these studies, in the discussions 
that Follow, we shall assume that the disk is not a con¬ 
tributing factor to the observed magnetic signatures. 


TABLE 1 

HD 30495 Properties 


Property 

Value 

Reference 

Spectral type 

G1.5 V 

(i) 

V 

5.49 ±0.01 

(2) 

B-V 

0.632 ±0.006 

(2) 

Parallax 

75.32 ± 0.36 mas 

(2) 

vsinz 

4.1 ± 0.8 km s- 1 

(3) 

T e fi 

5826 ± 48 K 

(4) 

log9 

4.54 ± 0.012 dex 

(4) 

[Fe/H] 

±0.005 ± 0.029 solar 

(4) 

Mass 

1.02 ± 0.01 Mq 

(4) 

Radius 

0.898 ± 0.013 Rq 

(4*) 

Luminosity 

0.837 ± 0.037 Lq 

(4*) 

Age 

970 ± 120 Myr 

(5*) 

Plot 

11.36 ±0.17 days 


A P 

0.59 ± 0.05 days 


An 

> 1.67 ± 0.15 deg/day 


sin i 

1.0 ±0.2 


i 

> 55.4° 


P eye, long 

12.2 ± [3.0] yr 


-^-cyc,long 

0.118 ± [0.044] 


-^short 

1.67 ± [0.35] yr 


-^short 

0.066 ± [0.028] 


Note. — 

References: (1") Gray et al. (20061 

(2) Hippar- 

cos, Perryman et al. (1997); van Leeuwen (2007) (3) Gai- 


"From BaumamTjgt al.| ( |2U10p measurements (5*) Derived 


from [Barnes| (|2007p age-rotation relationship. The lower 
section of the table are measurements found in this work 
Quantities in brackets represent one half of the observed 
range of values. 


Baliunas et al. 
variability Tn 


(1995) previously searched for cyclic 
30495 using the Lomb-Scargle pe- 
Mount Wilson S time series from 
Var”, defined as “sig- 
and 


HD 

riodogram on the 
1966-1992, but classified it as 
nificant variability without pronounced periodicity 1 
(as/S) > 2%. Tbe mean activity level was found to be 
high, ( S} = 0.297, as expected for fast rotators (compare 
to the so lar value from the same study ( Sq ) = 0.179). 
Hall et al. (2007) S'-index measurements from the Solar 
Stellar Spectrograph (SSS) again fou nd stronger-than- 
solar activity levels, with ( S) = 0.309. Hall et al. (2009) 
used twelve years of precise Stromgren b and y photom¬ 
etry from the Fairborn A utomated Photometric Tele¬ 
scope (APT) program (Henry 1999), finding a photo¬ 
metric variability for HD 30495 roughly six times solar. 
Furthermore, brightness was shown to decrease with in¬ 
creased chromospheric activity, indicating that variabil¬ 
ity in the star’s brightness is dominated by dark spots, 
typical of fast rotators in that study. 

Naively, and based solely on its similarity to the Sun 
and faster rotation and, hence, presumably greater dif¬ 
ferential rotation, we might expect HD 30495 to have a 
more vigorous dynamo, leading to higher magnetic ac¬ 
tivity and a shorter activity cycle. As we shall see, the 
former is borne out by observations, while the latter is 
not. 

2. ACTIVITY ANALYSIS 

We analyze a combined 47-year time-series of the 
Mount Wilson S'-index shown in Fig [Tj This dimen¬ 
sionless index is defined as the ratio of the core emis¬ 
sion in the Fraunhofer H and K lines of Ca II with the 
nearby continuum regions, as measured by the HKP-1 
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Fig. 1.— (a) HD 30495 combined S'-index time series, including data from MWO (•), SSS (o), SMARTS (■), CPS (a), and HARPS (▼), 
along with seasonal means (•), with error bars representing the error of the mean. The red curve is a 3-component sine wave model of the 
stellar cycles, while the horizontal red line is our reference point for global activity minimum, (b) Zoomed portion highlighting higher-cadence 
SMARTS data, with a P = 1.58 yr sine wave plotted for comparison (blue dashed curve), (c) APT differential photometry brightness 
measurements in the combined Stromgren b and y bands, in milli-magnitudes. Differences shown are HD 30495 nightly measurements (•) 
and seasonal means (•) with respect to the comparison stars, as well as the difference between the two comparison stars (□) (d) Seasonal 
mean differential brightness difference in the b and y bands, in milli-magnitudes, with colored regions indicating brighter b (blue) and y 
(green) emission. A horizontal dotted line indicates the grand mean in all panels. Magnitude scales are inverted such that brightness 
increases in the upward direction. The observations shown in this figure are available in the electronic version of this publication. See 
tables [3] and [1] in the Appendix. 
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and later HKP- 2 photometer s at Mount Wilson Q bser- 
vatory (MWO) (|Wilson||1968| |Vaughan et al.||1978|). Ca 
II H & K global-scale emission reversals are a signature of 
departure from radiative equilibrium, a defining feature 
of a chromosphere, and must be due to magnetic non- 
thermal heating mechanisms. Due to the unsurpassed 
duration and breadth of the Mount Wilson survey, S, 
a measure on an instrumental scale , has become the de 
facto standard for measuring stellar magnetic activity, 
and subsequent surveys have calibrated to the Mount 
Wilson scale. S is dependent on stellar properties such 
as temperature, surface gravity, and composition, which 
precludes its use for directly comparing activity levels 
of a heterogeneous ensemble. As we focus our analysis 
on a s ingle star, conversio n to a corrected quantity (e.g., 


R'hk > Noyes et al. (1984)) is not necessary. 


The combined^ time series in Figure [I]) a) contains 
1285 measurements from five different instruments. The 
majority of the measurements (624 measurements from 
1967-2003) come from the origi nal Mount Wilson sur - 
vey, calibrated as described in Baliunas et al. (1995). 
We assumed a uniform measurement error of 3% of the 
mean fo r this time series, near the upper limit quoted by 
Wilson (1968). The next largest portion of the measure¬ 
ments are from the Solar-Stellar Spectrograph (SSS) at 
Lowell Observatory (364 measurements from 1993-2014) 
(Hall & Lockwood 1995), taking the time series from the 
beginning of the Mount Wilson survey to present day. 
SSS obtains R ~ 10, 000 at Ha spectra, and S is de¬ 
rived by apeproximating the Ca II H, K and continuum 
bandpasses used by the MWO instrument. These data 
are then calibrated to the Mount Wilson instrumental 
scale using long-duration means of common targets. A 
typical measurement error of 2.4% was estimated using 
photon statistics in the K line core and detector prop¬ 
erties. Observations from the SMARTS Southern HK 
survey using the RC Spec R « 2500 spectrograph at 
1.5-m telescope at CTIO are the third largest contribu¬ 
tion (140 measurements from 2008-2013), and though of 
shorter duration, this queue-scheduled time series is not 
plagued by the large seasonal gaps of the other surveys, 
allowing short-period variation to be better determined. 
These data were calibrated to the Mount Wilson scale via 
common o bservations with SSS targets, as described in 
Metcalfe et al. (2010). An additional 108 measurements 
from 2011-2015 derived from HARP S R ss 120,000 spec¬ 
tra from a solar tw in planet search (|Ramfrez et al.||2014 
Bedell et al. 20151, again calibrated t o the MWO scal e 


using common targets, as described in (|Lovis et alTj~20 11 [) . 
Finally, we add 49 observations from 2002-2008 derived 
from R ss 55, 000 spectra of the Hamilton Spectrome¬ 
ter at Lick Observatory. These observation^] are part 
of the California Planet Search (CPS) and were simi- 
larly cali b rated to the MW O scale using common targets 
(|Isaacson & Fischer|[2010). 

Though each of these time series used a global calibra¬ 
tion to the Mount Wilson scale using long-term means 
of commonly observed targets, visual inspection of the 
combined time series revealed obvious discontinuities and 
differences in scale. This is likely due to the fact that the 
global calibration involves a compromise linear fit among 
all targets, while scatter about that fit reveals error in 

7 http://cdsarc.u-strasbg.fr/viz-bin/Cat?J/ApJ/725/875 


the calibration that would result in a discontinuity in any 
individual target. We applied a simple calibration that 
assumes overlapping periods of two different time series 
ought to agree on the mean for that period. To calibrate 
S to the scale of .So, the mean value over the period of 
overlap, S and Sq were calculated, and a scaling factor 
C = Sq/S was derived. The resulting calibrated time 
series S 1 = CS then has an equivalent mean value over 
the overlapping period to the base series Sq- The re¬ 
sulting scaling factors were C{ SSS —> MWO) = 1.015, 
C(SMARTS SSS') = 1.067, C(HARPS SSS') = 
1.074, C(CPS -A SSS') = 1.098. The SMARTS, HARPS, 
and CPS time series were scaled using overlapping por¬ 
tions of the post-calibration SSS time series, therefore 
their overall scaling is multiplied by C(SSS -A MWO). 
This calibration removed obvious discontinuities in the 
combined time series and reduced the standard devia¬ 
tion by 3.8%. The final combined time series has a grand 
mean S = 0.303 and a standard deviation a = 0.0167. 
Seasonal means for the combined time series are shown 
as black circles in Figure [TJa). Following these seasonal 
means, clear cyclic behavior is visible, emphasized by the 
cycle model (red curve) described below. 

We also examined the 22-yr time series of differential 
photometry acquired with the T4 0.75 meter Automatic 
Ph otoelectric Telescope (APT) at Fairborn Observatory 
(Henry 1999), shown in Figure flj(c). These measure¬ 
ments, made in the Stromgren b (467 nm) and y (547 
nm) bands, are a difference with respect to the mean 
brightness of two stable comparison stars, HD 31414 and 
HD 30606. The differential measurements in the b and y 
bands are then averaged to (b+y)/ 2 to create a “fry” band 
that increases the signal to noise ratio. The unimpor¬ 
tant mean difference is subtracted from the time series. 
The stability of the comparison stars is demonstrated in 
the seasonal mean of their brightness difference in the by 
band, shown as white squares in FigurefHc), with a stan¬ 
dard deviation cr = 0.00093 mag. HD 30495 by brightness 
is strongly variable (cr = 0.0065 mag) and out-of-phase 
with the S'-index shown in Figure [It b). A rank correla¬ 
tion test between S and by seasonaFmeans shows 99.98% 
significance in the correlation, which is plotted in Figure 
[2](a). This is interpreted as evidence the star’s brightness 
variations are dominated by dark spots, which are more 
prevalent during times of activity maximum. Figure[l](d) 
plots the A(b - y) color index where blue shading indi¬ 
cates negative color index and green shading indicates 
positive color index. Comparing panels (b), (c) and (d) 
of Figure [T| we see that HD 30495 gets bluer as it gets 
brighter (activity minimum) and redder as it gets fainter 
(activity maximum). This is shown again more clearly in 
Figure [2] (b) and (c), in particular the remarkably tight 
color-brightness correlation. We interpret this color shift 
as an increase in surface temperature during times of ac¬ 
tivity minima, due to the reduction of cool spots on the 
surface. 


We c omputed the Lomb-Scarg le periodogram (Scargle 
(1982), Horne & Baliimas| (| 1 986|) from our time series to 
find statistically significant periodicities in the data, with 
the results shown in Figure [3] To verify the robustness of 
the peaks, we compare the periodogram of the combined 
time series (thick black line) to those of the individual 
MWO, SSS, and SMARTS series over shorter intervals, 
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Fig. 2.— Correlations among seasonal means of activity, brightness and color from the time series of Figure [T] Error bars indicate the 
error in the mean. Magnitude scales are in milli-magnitudes and are inverted such that brightness increases in the upward/rightward 
directions. 


as well as to the periodogram of the (b+y )/2 photometry. 
Not shown in the figure are the large peaks beyond 25 
years in the MWO and combined S'-index periodograms, 
which are most likely due to the windowing of the entire 
time series, not true physical variation. The hatched re¬ 
gions of P < 1.1 on the left side of periodograms contain 
a number of large peaks near 1 year, which are aliases 
due to the seasonal sampling in our time series. We ver¬ 
ified these are all aliases by obtaining a least-squares fit 
of the data to a sine wave with a period set by one of the 
~1 yr peaks, then subtracting that signal from the time 
series and re-computing the periodogram. The new pe¬ 
riodogram would no longer contain the ~1 yr peak, and 
a corresponding low-frequency peak would be removed 
as well. This established a symmetry between the low- 
frequency peaks and these ~1 year peaks, and as a result 
we do not consider any peak <1.1 years to be physical. 
(See also Figure^c), in which spurious ~1 yr peaks are 
found in the periodogram of a signal of pure sine waves 

of lower freq uency.) _ 

Following Horne & Baliunas fl986|, we calculate the 
“False Alarm Probability" (FAF) threshold: 


z = -ln(l- (l-F) 1/iVi ) 


( 1 ) 


where F is the probability that there exists a peak 
of height z at any frequency due to random Gaussian 
noise in the signal, and iV, is the number of independent 
frequencies in the time series. We computed A/ by gen¬ 
erating 5000 random time series with the same sampling 
times of our data, generating a probability distribution 
for the maximum peak z , and fitting this distribution to 
equation 0 inverted for F, with A/ as the free parame¬ 
ter. The upper threshold (green line) shown in Figure [3] 
corresponds to F = 10~ n , the th reshold for an “excel¬ 
lent” cycle in |Baliunas et al. (1995), and the lower thresh¬ 
old (red line) is for F = 10“^ (99.9% significance), the 
minimum requirement for a “poor” cycle in that work. 
The FAP thresholds shown are those computed for the 
combined time series, however they are similar to those 
obtained for the individual component time series, being 
only slightly more stringent. 

The uncertainties in peak positions were estimated us¬ 


ing a Monte Carlo method. In each trial, each time series 
measurement is randomly sampled from a Gaussian dis¬ 
tribution defined by that measurement’s value and un¬ 
certainty. Then, a periodogram is computed and the new 
peak position saved. By running 5000 trials, an approx¬ 
imately Gaussian distribution of peak positions is ob¬ 
tained, and the uncertainty is estimated as its standard 
deviation. 

In the combined time series we found four resolved 
peaks above the “excellent” threshold: a long period 
peak Piong = 12.77 ± 0.09 yr, and a cluster of three 
short-period peaks at P s hort,i = 1.572, P s hort,2 = 1-486, 
and P s hort,3 = 1-615 yr (cr sho rt ~ 0.003). The Pi on g peak 
is found between nearby significant peaks found in the 
MWO (Piong, MWO = 10.7 yr) and SSS (P long ,sss = 15.3 
yr) time series, and nearby the Pi on g,APT = 12-2 peak 
from by photometry. The spread in periods from the ear¬ 
lier MWO data to the later SSS data indicates that, like 
the Sun, the long-term cycle is only quasi-periodic, and 
the duration of each individual cycle varies. The ~13 yr 
peak has a protruding “shoulder” on its right side, which 
is due to an unresolved peak near 17 yr. This peak was 
resolved in ~25% of the Monte Carlo trials done to de¬ 
termine the uncertainty in peak positions, allowing us to 
measure a mean value of 16.9 yr. This ~17 year period 
and another large peak at ~37 yr were found to be ar¬ 
tifacts of the amplitude structure and/or the duration 
of the time series (data window). We verified this by 
computing periodograms of various fractions of the data 
window (e.g. 2/3 to 1/2 of the total duration) at various 
offsets and noting that the ~17 yr peak disappears in all 
cases and the ~37 yr peak shifts close to the duration of 
the new window. 

Figure [3|b) shows the short-period peaks are almost 
perfectly matched by two peaks (1.49 yr and 1.61 yr) in 
the APT periodogram. We also find a broad correspond¬ 
ing “excellent”-class peak at 1.63 yr in the SMARTS time 
series and again in the HARPS series at 1.75 yr. A less 
significant “poor”-class peak at 1.85 yr is in the SSS data 
and 1.53 yr in the MWO data. The spread in values in¬ 
dicates that these short-period variations are not of a 
constant frequency, which we investigate in detail in Sec¬ 
tion [4] 
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Period [yr] 


Fig. 3. — Lomb-Scargle periodograms from the time series of Figure [I] Panel (a) contains the result from single-instrument S-index 
surveys and panel (b) the combined S-index time series, as well as the APT photometry of Figure me) . Note the division in the period 
scale. The hatched region near P « 1 yr contains artifacts of the seasonal sampling. The green ana red ho rizon tal dashed lines are the 
“excellent” and “poor” significance thresholds for the S-index periodograms, as defined in Baliunas et ah (| ] 995|i. Note that the APT 
periodogram is scaled down by a factor of five for easy comparison and the magenta horizontal line is the “excellent” threshold for that 
series. 


We find that both the long and short-period signals are 
found consistently in several distinct .S-index time series 
of different time intervals, as well as in the APT differen¬ 
tial photometry, a measurement using very different ob¬ 
servation methods to sample a physically distinct region 
of the star. This is strong evidence for the co-existence of 
variability on different time scales and in distinct regions 
of the stellar atmosphere, analogous to the solar observa¬ 
tions of the 11-year cycle and quasi-biennial oscillation. 

3. SIMPLE CYCLE MODEL 

In the Sun, each occurrence of rising and falling activ¬ 
ity is numbered and, somewhat confusingly, referred to 
as a “cycle”, with the current episode denoted as “cycle 
24”. Properties of each cycle such as duration, ampli¬ 
tude, and shape are measured and found to vary. We 
wish to similarly decompose the ~12 yr periodic signal 
of HD 30495 into individual cycles and measure their 
properties. For the Sun, this decomposition is typically 
done by identifying cycle minima in a smoothed time 
series of a proxy such as sunspot number, e.g. the 13- 


month boxcar smoothing of the monthly average s, and 


then using the minima as delimiters for each cycle (Hath¬ 


away 


2010). The seasonal gaps in stellar time series do 


not allow us to use this same prescription. Instead, we 
construct an idealized smoothed model of the time series 
as a superposition of low-frequency sine waves: 


S(t)i = A, sin (t + + y t (2) 

where Pi is a low-frequency period from a periodogram 
analysis, and the amplitude A,, phase <f>i and offset yi are 
found using a least squares optimization of this model to 
the mean-subtracted data. The final model is simply: 


N P 

S(t) = J2S(t)i + S (3) 

i 

where Np is the number of component sine waves and 
S is the original mean S-index. 
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Period [yr] 


Fig. 4.— Lomb-Scargle periodogram of (a) the original composite 
iS-index time series (b) 3-component cycle model with equal-spaced 
sampling (c) 3-component cycle model with the same sampling as 
the original data (d) the residual of the original data minus the 
cycle model. Note that period and power scales change at the 4-yr 
mark; the left and right y-axis give the power scale for that side. 
Periodograms of the cycle model are normalized to 1. 


TABLE 2 

HD 30495 Cycle Properties 


Cycle 

Start 

Max 

Duration 

Smax 

Acyc 

0 

(1961.7) 

1969.0 

(15.3) 

0.324 

0.156 

1 

1977.1 

1981.2 

9.6 

0.297 

0.067 

2 

1986.7 

1993.7 

11.7 

0.305 

0.095 

3 

1998.4 

2005.8 

15.5 

0.324 

0.156 

4 

2013.9 

- 

- 

- 

- 


We obtain the parameters of equation i> by iteratively 
finding the lowest-frequency period in a Lomb-Scargle 
periodogram of the composite S'-index time series, fitting 
the sine to the data using a least-squares optimization, 
and subtracting the result from the data before comput¬ 
ing the next periodogram. We carry out three such itera¬ 
tions in order to find the three low-frequency components 
revealed by the periodogram of the original dataset. The 
resulting model is akin to the result of a low-pass filter on 
the data. Indeed, subtracting the model from the data 
and computing the periodogram reveals that we have ef¬ 
fectively removed most of the low-frequency power, as 
can be seen in Figure [4] The parameter set obtained 
from the fits of equation |2| were P = {37.6,18.8,12.2}, 
A = {11.4,4.68,6.64} x 1(T 3 , <j> = {—0.103,5.93,5.45}, 
and y = {-2.47, -0.0519,1.26} x 10" 3 . 

Using this model we characterized the individual cycles 
of HD 30495 as shown in Table [2j Quantities in paren¬ 
thesis are based on an extrapolation beyond the data, 
and should be treated with caution. We see a spread in 
cycle durations from 9.6 to 15.5 yr, with a mean value 
of the fully observed cycles 1-3 of 12.2 yr. This value is 
close to the Pi ong = 12.77 yr peak in the periodogram of 


Figure [3jb) . We shall adopt the mean value of 12.2 yr 
as our Rest estimate for the mean cycle period for HD 
30495. The spread of individual cycle durations from 
the model gives us an estimate of the variability of cycle 
durations, APi ong /2 = 3.0 yr. The increase in cycle du¬ 
ration from cycle 1 to 3 obtained from our model agrees 
with the trend observed in the periodograms, with the 
earlier MWO series which includes cycles 1 and 2 having 
Pong, mwo = 10.7, and the later SSS series which only 
includes cycle 3 having Pi ong ,sss = 15-3 (see Figure[3](a)). 
The cycle model components at P = 37.6 and 18.8 yr are 
artifacts of the amplitude structure in the time series and 
the data window, as discussed above. These components 
are necessary, however, to reconstruct the amplitude and 
cycle duration variability present in the original time se¬ 
ries. For example, the amplitude A = 6.64 x 10 -3 from 
the P = 12.2 yr least-squares fit to the original data is 
not representative of the amplitude of the individual cy¬ 
cles seen in Figure [lja), being at least a factor 2 too low. 
Only when the other two low-frequency components are 
included does the model amplitude more closely match 
the data. 

The relative cycle amplitudes from the model are 
shown in column 6 of Table [ 2 J defined as A rvr = A S/S 
follow ing |Soon et al.| (|1994D and |Saar & Brandenburg 
(20021 studies of cycle amplitudes from the Mount Wil- 


son program. Determining AS = S max — S le { requires us 
to select a reference point that approximates the lowest 
possible activity level for this star. Our cycle model does 
not effectively model the depth of the minima, as can be 
seen in Figure [l]( a), which precludes its use for setting 
S le f. Instead, we choose S re { to be the lowest seasonal 
median with more than 3 measurements, S le { = 0.2763 
from the 1986-1987 season. This definition of activity 
minimum avoids too much sensitivity to outliers, while 
being low enough that only 3.6% of the data lie below 
this point. With our choice of S re f and S = 0.303 we 
find A cyc ranging from 0.067 to 0.156, and a mean am¬ 
plitude A cyc = 0.118. We find an increase in cycle am¬ 
plitude from cycles 1 to 3, which occurs in parallel with 
the increase in cycle duration. Note that for the Sun the 
amplitude of a subseq uent cy cle is negatively correlated 
with cycle period (Hathaway (2010). The two transitions 
of fully observed cycles tor HI) 30495 indicate a positive 
correlation, though obviously no firm conclusions can be 
drawn from so little data. 

Robustness of the cycle model was examined by run¬ 
ning Monte Carlo simulations with the nightly measure¬ 
ments resampled from within their estimated uncertain¬ 
ties for each trial, then repeating the iterative procedure 
described above to compute the cycle model periods, am¬ 
plitudes, phases, and offsets. The resulting minima posi¬ 
tions, maxima positions, and amplitudes were gathered 
and the standard deviation of those distributions com¬ 
puted. The times of minima and maxima were found 
to be relatively robust forming a Gaussian distribution 
with a standard deviation of a = 0.25 yr. The amplitudes 
were also fairly robust, forming a Gaussian distribution 
with er = 8.6 x 10 -4 . We use this to estimate the un¬ 
certainty in AS to be about 7%, which dominates the 
uncertainty of our relative amplitudes A cyc = A S/S. 


4. SHORT-PERIOD TIME-FREQUENCY ANALYSIS 
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Fig. 5.— Short-Time Lomb-Scargle analysis. A Lomb-Scargle periodogram is computed for 5-yr moving windows in 1 yr increments of 
the composite 5-index time series for HD 30495. The contour plot gives the periodogram power normalized by the number of data points 
as a function of time and period, with the 99.9% significance contour highlighted as a dotted white line. The indicated time is for the 
center of the 5-yr window. In each window, the highest peak above the 99.9% threshold is found and plotted as a large open circle, while a 
secondary significant peak, if present, is plotted as a small open circle. Vertical blue and red dashed lines indicate the minima and maxima 
in the long-term cycle model.The top panel indicates the number of data points in each 5-yr window. The right panel gives the integrated 
normalized power for all windows. The bottom panel plots the amplitude A5 of sine fits of the significant peak periods in the windowed 
data as open circles, with the amplitude of the 3-component cycle model in red. 


A sine fit with P = 1.58 was found for the combined 
time series for a 5-year window centered at 2010.0 is 
shown in Figurejljb), with the curve colored cyan within 
the fitting window and colored gray outside of the win¬ 
dow. Comparing this sine wave to the data visually 
demonstrates the short-period variation, especially vis¬ 
ible during the epoch of SMARTS measurements (cyan 
points). Even in this short segment we notice that the 
data goes out-of-step with the sine curve after 2012. This 
shows by example the quasi-periodic nature of the short- 
period variations, which was also evident in the triplet 
of peaks near P « 1.6 yr in Figure[3](b). 

To investigate possible coupling between the short- and 
long-period variations found above, we performed a peri¬ 
odogram analysis on a 5-yr moving window of our com¬ 
bined S'-index time series, the results of which are shown 
in Figure [5j First, we removed the low-frequency com¬ 
ponents from the data by subtracting the cycle model 
described in Section [3] Comparing Figure J4] panels (a) 
and (d), we see that this procedure amplifies the high- 
frequency peaks in the residual periodogram. However, 
we found the results of this analysis are qualitatively the 
same when working with the original composite time se¬ 
ries. Next, the residual time series is divided into 1-yr 
bins, which contain the entirety of the MWO/SSS sea¬ 
sonal observations. Every set of five consecutive bins are 
then subjected to a Lomb-Scargle periodogram analy¬ 
sis as described in Section [2j searching for periods from 
1.1 to 2.5 years. Five year windows were chosen as a 
sufficient time period to capture multiple oscillations of 
the P ss 1.6 yr variability detected in Figure [3j)b) . The 
period search cutoff at 1.1 yr is to avoid aliasing issues 


related to the seasonal sampling, and the 2.5 yr cutoff 
is to avoid strong signals associated with the duration 
of the window. In total, 44 periodograms were calcu¬ 
lated as well as a 0.001 (99.9% significance) false alarm 
probability threshold based on a Monte Carlo analysis as 
described in Section [2] Up to two significant peaks were 
extracted from these periodograms, and a least-squares 
fit of a sine function to the data window is done to de¬ 
termine the amplitude of that signal. In Figure [5] the 
normalized periodogram power is plotted as a contour 
area plot in the time-period plane, and the positions of 
significant peaks are plotted as open data points. The 
99.9% significance contour is plotted as a dotted white 
line. 

It is important to note that higher periodogram power 
is possible in time series with a larger number of data 
points, so that periodogram power does not scale equally 
from one window to the next. To check this behav¬ 
ior, we generated unevenly sampled time series of a sine 
function by randomly removing 50% of the data points 
and calculated the Lomb-Scargle periodogram, finding 
that periodogram power is proportional to the number of 
data points. To correct for this, we normalized the peri¬ 
odogram power by the number of data points in the win¬ 
dow. Differences in periodogram power between windows 
are then due to differences in amplitude of the underly¬ 
ing signal. The number of data points in each window is 
plotted in the top panel of Figure [5] for reference. The 
99.9% significance threshold is computed separately for 
each window, so the contour position is correct despite 
the difference in the number of data points. 

Figure [5] reveals that variability near P s hort ~ 1-6 yr is 
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intermittent with peak periods occuring across the full 
range of periods analyzed. Significant peak variability in 
the 1.4-1.8 yr range begins in 5-yr windows centered at 
1985, 1993, 2000, and 2009, with each episode lasting for 
3-5 years. The 1985 and 2009 episodes precede the cy¬ 
cle 2 and 4 minima in the long-period cycle, respectively, 
while the 1993 and 2002 episodes roughly coincide with 
cycle maxima. In contrast to the minima of cycles 2 and 
4, the cycle 1 and 3 minima are devoid of short period 
variability. Peaks are often found from ~2.2 to ~2.4 yr 
as well, which is confirmed in the plot of the power inte¬ 
grated over all windows. Here the broad peak at ~1.7 yr 
is notably shifted from the cluster of narrow peaks in the 
periodogram of Figure [3jb) . This may be due to vary¬ 
ing phase in intermittent short- period signals leading to 
interference effects in the full-time-series periodogram. 
We will take the mean of the detected peak periods < 
2.0 yr as our best estimate of the short-period signal, 
Pshort = 1-67 ± [0.35], where the quantity in brackets is 
half the observed range of peak periods. 

We analyzed the peak-to-peak amplitudes AS of sine 
fits to the data with the significant peak periods in the 
range of 1.1 to 2.0 yr and found them to range from 0.012 
to 0.030, with a mean of 0.020. The average short-period 
relative amplitude A s hort = A S/S is then 0.066 ± [0.028], 
roughly half of the average long-period amplitude but 
nearly equal to relatively low amplitude of cycle 1, as 
deduced from our cycle model. We performed a rank- 
correlation test between the short-period amplitudes and 
the long-period cycle model, but no significant correla¬ 
tion was found. 

From the above observations we conclude that there is 
no clear association between the long-period cycle and 
the episodic short-period variations. The presence or 
absence of the short-period variations are found in all 
phases of the long-term cycle, and the amplitudes are 
not correlated with the long-term cycle amplitude. 

5. ROTATION 


in|Gaidos et al. (2000) using 
I. The dense sampling of the 


We repeated the rotation measu rements done for six 
seasons of APT photometry : 
the current 22-season record. 

APT program allows the detection of rotational modu¬ 
lations due to the transit of spots on the stellar photo¬ 
sphere. The time series is broken into individual seasons 
containing 55-185 measurements over the course of 150- 
200 days. From each season’s time series a Lomb-Scargle 
periodogram was computed looking for rotational peri¬ 
ods from 2-25 days. Peaks passing the 99.9% significance 
level are taken as a signal of rotation. Uncertainty in 
the peak position was calculated using the Monte Carlo 
method described in Section [2] and ranged from 0.014 to 
0.2 days. In total, rotational periods were detected in 
17 of 22 seasons. Taking the mean, we find the average 
Prot = 11.36 ± 0.17 d. This resu lt is within the uncer¬ 
tainty of the previ ous Gaidos et al. d2000 |) measurement, 
as well as that of iBaliunas et ai.l (19961 who obtained 


rotation from a densiy-sampied season of MWO S'-index 
observations. The range of the precisely determined ro¬ 
tations are from 10.970 ± 0.028 to 11.560 ± 0.023 days, 
giving A P = 0.59 ± 0.05 d. This is reduced from the 


Gaidos et al. (2000) value of A P = 1.0 d, due to the fact 
that we could not reproduce a significant 10.5 d detec¬ 
tion in season 6. AP is interpreted as a sign of surface 



Fig. 6. — Top two panels: time series of seasonal rotation pe¬ 
riod measurements beneath seasonal mean S-index time series, for 
comparison. Filled (open) circles are measurements for seasons in 
which a rotation period detection was successful (unsuccessful). A 
horizontal dotted line marks the grand mean for the whole time 
series. Bottom panel: seasonal activity-rotation correlation plot, 
with error bars representing 1 -<j uncertainty in the rotation period 
and seasonal mean S value. Data points are annotated with the 
two-digit year. 


differential rotation, due to transiting spots at various 
latitudes. Our measurements range over one and a half 
stellar cycles, which increases the confidence that we have 
explored the full range of rotational periods which can 
be sampled with this technique. However, due to the 
unknown latitude ranges of the spots on the star, the 
measurement provides only a lower bound on the true 
equator-to-pole surface differential rotation. 

The time series of rotation measurements and the sea¬ 
sonal activity-rotation relationship are shown in Figure 

t P ro t vs. S' is a kind of pseudo-butterfly diagram in 
e absence of spot latitude information, which would 
demonstrate diffe rent morphologies under different mi¬ 
gration patterns (Donahue 1993). For example, in the 
Sun, if rotation were measured by tracking spots during 
the cycle we would expect a long rotation period and low 
activity at solar minimum (high latitude spots), transi¬ 
tioning to shorter rotation periods and higher activity 
until solar maximum (mid-latitude spots), and finishing 
with still shorter rotation periods as activity wanes (near- 
equa tor spot s). A variety of morphologies were observed 


Donahue (1993), including the anti-solar case. From 


Figure [6] we see that for HD 30495 in general rotation is 
slow when activity is high (with the exception of the 2008 
season), but when activity is low both long and short ro¬ 
tation periods are seen. Tracing points in chronological 
order reveals no clear pattern, with seasons transition¬ 
ing from quadrant to quadrant in the diagram. Near the 
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maxima of cycle 3 (t = 2005.8, ref Tableful) there is a clus¬ 
ter of four seasons with slow rotation, 2004-2007. The 
1997-1998 rotation was also slow, which occurred just 
before the minima at the start of cycle 3 (t = 1998.4), 
however the 2011 and 2013 rotation before the start of 
cycle 4 (t = 2013.9) was relatively fast. The lack of co¬ 
herent structure in Figure [ 6 ] may be an indication that 
large spots which make the detection of P rot possible are 
not restricted to a narrow range of latitudes as for the 
Sun. 

Line-of-sight inclination is an important factor for in¬ 
terpreting the stellar rotation and activity data. We es- 
timated the i nclina tion using the v sin i measurement of 
Gaidos et al. (2000), together with v = 2nR/P Iot , where 
R was derived using a oc M/R 2 and th e mass and surface 


gravity estimates of Baumann et al. (2010) (see Table 
fy, as well as our rotation measurement! We obtained 
smi = 1.0 ± 0 . 2 , indicating an equator-on view of the 
star, but the large uncertainty giving a one-sigma range 
of i > 55.4°. This perspective provides a best-case sce¬ 
nario for measuring rotation with spot transits, as well 
as a “solar-like” view of the activity cycles. 

Finally, using t he age-rot a t ion relationship given in 
equation (3) of Barnes (2007) along with our mean ro¬ 
tation period, we obtain an age of 970 ± 120 Myr. This 
revises the f gyr0 from Table 3 of that work, which was 
based on a lower estimate for the rotation period. 

6 . DISCUSSION 

We have observed quasi-periodic signals with repre¬ 
sentative values of Pi ong = 12.2 ± [3.0] yr and P s hort = 
1.67 ± [0.35] yr in the chromospheric activity of the fast- 
rotating (P rot = 11.36 ± 0.17 days) solar analog HD 
30495. A simple three-component sine cycle model shows 
three full cycles in the time series, each with varying du¬ 
ration and amplitude. This, combined with the improved 
signal-to-noise from our longer time series, has allowed 
us to demonstrate the “pronounced periodicity” that was 
previously lack i ng to classify this star as cycling in the 


Baliunas et al. (1995) study. 

Taking the ratio of these periodicities to the rotation 
period, n = P C y C /Prot, we find ni ong ~ 400 rot/cyc and 
Jtshort ~ 50 rot/cyc respectively, which closely corre¬ 
spond to the “active ” and “inactive” seque nces found in 
the stellar sample of Bohm-Vitense (2007). What is re¬ 
markable from that work is that the Sun appeared to be 
a unique outlier in that sample, with ni ong0 ~ 150, and 
taking the solar QBO period as 2 years, n s hort ,0 ~ 30. 
Why does the Sun appear as an outlier with respect to 
the stars? One explanation could b e that the relatively 
small sample of 21 stars in Bohm-Vitense (2007) is in¬ 
sufficient to show the full picture of the relation between 
P cyc and P ro t) and further data will simply erase the 
observed trend s. Indeed, the sample of Sa ar fe Bran-| 
denburg (1999), which is a superset of the Bohm-Vitense 
sample, includes a few neighboring points for the Sun 
on log plots of quantities proportional to P cyc and P ro t- 
Our datum for HD 30495, however, is decidedly not 
a neighbor of the Sun on these plots, while it agrees 
well with the trend set by stars on the active branch. 
Both the Sun and this star have a similar time scale for 
the observed long and short periodicities, yet rotation, 
n*/f2 0 « 2.3, is very different. This poses a serious 


problem for Babcock-Leighton flux transport dynamo 
models whose time scale is determined by the meridional 
flow. 3D hydrodynamic models show that meridional 
flow speed decreases with increased rotation rate (Ballot 


et al.|2007 |Brown et al.|2008 ), but kinematic mean-held 


dynamos with one meridional flow cell have cycle times 
proportional to the fl ow speed a nd hence slower cycles 
with faster rotation. Jouve et al. (2010) investigated this 
problem, finding that multiple? meridional flow cells in 
latitude were needed to make cycle period proportional 
to rotation. Unfortunately, observations to determine 
whether this indeed occurs in fast-rotating stars such as 
HD 30495 are nearly beyond imagination due to the ex¬ 
tremely slow flow speeds, of order ~10 m/s on the Sun. 

We also compared relative amplitudes A cyc = A S/S of 
t he lo ng an d short variatio ns to the stellar ensemble in 


Saar & Brandenburg (2002). We found the average cycle 
amplitude to be lower than the lowest star in that sample, 
though not excessively far from other active-branch stars, 
while the maximum amplitude (cycle 3) matches closely 
with the other active branch stars of similar B — V color. 
That the average val ue is low is likely du e to di fferences in 


methods used by Saar & Brandenburg 


et al. (1994), the origin of part of the sample) which 


| 2002 | ( and Soonj 


estimated peak-to-peak amplitude using the entire 25- 
year time series of the MWO sample. In addition, our 
method used the 3-component cycle model and a global 
minimum defined by the lowest seasonal median in the 
time series, which by design filters out the variance from 
high-frequency components, resulting in a lower estimate 
in amplitude. 

HD 30495 is nearly rotation and cycle-degenerate with 
the K2 star e Eridani, which has P ro t = 11.1 day s and 


Pinny = 12 .7 yr, but a longer P s hort = 2.95 yr (Met¬ 


calfe et al. 20131. In this comparison, the experimental 


variables include the depth of the convection zone, with 
stellar structure models predicting a deeper convective 
region for the cooler K star, as well as the average con¬ 
vective velocity. Apparently, these factors alone are not 
enough to prescribe a substantially different long-term 
cycle period. 

In most dynamo theories, differential rotation is the 
driving force of the D-effect, responsible for turning 


poloidal magnetic flux into toroidal flux (e.g. Babcock 


1961). Stars with greater differential rotation would be 


expected to “wind up” the field faster, and we should rea¬ 
sonably expect a shorter cycle time, as at least half of the 
process would be faster (the other half being the return 
of toroidal field to poloidal). Using our range of rotation 
period detections, the total measured surface rotational 
shear is AD = 1.67 ± 0.15 degrees/day. Then , using the 
the solar s urface differential rotation result of|Snodgrass| 


& Ulrich (1990), we calculate the equator-to-poie total 
shear for comparison, finding AD*/AD 0 > 0.40 ± 0.03, 
which is a lower bound due to the unknown latitude 
ranges causing the rotation signal on the star. However, 
if one is prepared to assume that the spot latitudes of 
HD 30495 never form above 45°, as for the Sun, then 
differential rotation can be compared in terms of the so¬ 
lar shear from the equator to 45°, giving Af 2 */AD, 45 i0 > 
1.02 ± 0.09. Both results allow the possibility that llD 
30495 and the Sun have equivalent surface differential 
rotation, which might help to explain their similar cy- 
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cle characteristics. Asteroseismic measurements of HD 


30495 may be able to put tighte r limits on Afl (Gizon & 


Solanki 2004 Lund et al. 20141. 


(Jur time-frequency analysis does not indicate coupling 
between the intermittent short-period variability and the 
long-period cycle for this star. This is in contrast to the 
solar QBO, wh ose amplitude is strongly modulated by 
the 11-yr cycle ( Bazilevskaya et al.|2014 ). For HD30495, 
the amplitudes of the short-period variations are at times 
larger than the amplitude of the long-term cycle, bring¬ 
ing into question which component of the variability is 
more fundamental to its dynamo. The absence of cor¬ 
relation between the two periodicities may be an indi¬ 
cation that they are of a fundamentally different na¬ 
ture. It would be interesting to know at what time 
scale global magnetic field polarity reverses in this star, 
if indeed it does reverse. The S time series during the 
densely-sampled SMARTS era shows a convincing sinu¬ 
soidal variation, perhaps as convincing as the long-term 
trend. A campaign of Zeeman Doppler Imaging mea¬ 
surements spaced over at least a four year period should 
be able to determine if the short-period variability is 
polarity-reversing as well. 

To the best of our knowledge, this is the first work to 
separate and characterize the amplitudes and durations 
of individual cycles from a stellar activity proxy. There 
are rich opportunities in this direction to explore the vari¬ 
ability of the cycles themselves , as well as differences in 
stellar behavior during times of minima and maxima, 
which in turn can provide additional constraints for dy¬ 
namo models. Already, the periodic signals measured 
here, together with the global properties collected in Ta¬ 
ble [l] present a well-characterized object to study with 
dynamo models. The existence of two significant time 
scales of variability in activity poses an additional mod¬ 
eling constraint. This bright object is also a prime can¬ 
didate for future asteroseismic observations, which can 
further constrain its mass, rad ius, rotat ion profile, and 


depth of the convection zone (Metcalfe 2009). Success¬ 
ful modeling of such well-described targets will hopefully 
lead to improved understanding of the dynamo process. 
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APPENDIX: OBSERVATIONS 

The observational data used in this study are available 
in the online version of this publication. Table J3] shows 
a sample of the nightly S measurements available from 
MWO, SSS, SMARTS,' CPS, and HARPS shown in Fig¬ 
ure [I] (a) and (b). Note that the equal-mean calibration 
described in Section [2] is applied to these data; the data 
can be returned to its original calibration by dividing by 
the constants described there. Table [4] gives a sample of 
the nightly differential photometry measurements in the 
Stromgren b and y bands. Letting Ci be the itli column 
of Table|4j the A(6 + y)/2 time series of HD 30495 in 
Figure [2](cJ is obtained using: 

A(t + 9 )/2„ D3M8s = i(^ + 54^) 

and subtracting the mean. The difference between the 
comparison stars is given by: 

A (b + y)/ 2 comps = C& + C ‘ 

The seasonal means of this series are shown as white 
squares in Figure[l|c). Finally, the color difference series 
A (b — y) of Figurejljd) is given by: 

A /, \ C 2 + C 4 C 3 + C 5 

A(0 - 1/)hD30495 — -^- 9 - 
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